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We apply a general method for the estimation of com- 
pletely positive maps to the l-to-2 universal covariant cloning 
machine. The method is based on the maximum-likelihood 
principle, and makes use of random input states, along with 
random projective measurements on the output clones. The 
downhill simplex algorithm is applied for the maximisation of 
the likelihood functional. 



PACS Numbers: 03.65.-w, 03.67.-a, 03.67.Lx 



In this paper, we apply the general method of Ref. 
to the l-to-2 universal covariant cloning machine for 
spin- 1/2 systems. In Section II we briefly review the 
ML principle and its use in measuring quantum devices. 
Section III presents the calculations needed to explic- 
itly construct the likelihood functional to be maximised. 
Some numerical results obtained through a Monte Carlo 
simulation and the simplex searching algorithm are then 
shown to confirm the reliability of the method. Section 
IV is devoted to conclusions. 



II. MEASURING QUANTUM DEVICES 



> 

o 
o 
o 

o 
o 

^ : 

-i— > ■ 
G ■ 
cd 

s : 
cr 



I. INTRODUCTION 

Perfect cloning of unknown quantum systems is for- 
bidden by the laws of quantum mechanics Q. However, 
universal covariant cloning has been proposed 2 ] , and has 
been proved to be optimal in terms of fidelity 3 ^ . Some 
relevant applications of cloning is eavesdropping in quan- 
tum cryptography B, and engineering of new kinds of 
joint measurements |5| . Now suppose that someone pro- 
vides you with a physical device, telling you that it works 
as a universal cloning machine. How can you check ex- 
perimentally this statement? You could perform different 
kinds of measurement on the output of the device, along 
with different preparations of the input states. Then, you 
could compare the correlations you have inferred from 
the outcomes with your theoretical predictions. How- 
ever, the most reliable and general way to proceed could 
be entirely reconstructing the completely positive (CP) 
map that univocally characterizes a physical device. 

Recently, a general method to solve this relevant issue 
has been proposed in Ref. . The method is based on the 
maximum-likelihood (ML) principle, applied to the data 
obtained by random measurements on the output of the 
device. The ML method has been used in the context 
of phase measurement (gj, and to estimate the density 
matrix ||, some parameters of interest in quantum op- 
tics [ [Tof , and the CP maps of quantum communication 
channels Q. For the problem of CP map reconstruction 
, the constraints inherent the maximisation can be im- 
posed by exploiting the isomorphism between CP map 
from Hilbert spaces Ti. to K. and non-negative operators 
in the tensor-product space K.®H (11]-(l| • Such isomor- 
phism has already been useful for the study of positive 
maps Ej| and to address the problem of separability of 
CP maps El. 



The maximum-likelihood principle states that the best 
estimation of unknown parameters is given by the values 
that are most likely to produce the data one experimenter 
has observed. Hence, this principle involves the maximi- 
sation of a function of the unknown parameters that is 
given by the theoretical probability of getting the col- 
lected data. 

Consider a sequence of K independent measurements 
on the output of a physical device acting on quantum 
states. Each measurement is described by the element 
Fi(xi) of a POVM, where xi denotes the outcome at the 
Ith measurement, and I = 1,2, ...,K. Let us denote by 
pi the state at the input at the Ith run. The probability 
of getting the string of outcomes x = {x%,X2, ...,xk} is 
given by 



p(5f)=n£ =1 Tr[£(p l )F l {xi)] 



(1) 



The best estimate of the map £ maximizes the logarithm 
ofEq. ffih 



K 



C{£) =Y J ^g^[£{pi)F l {x l )\ 



(2) 



i=i 



over the set of completely positive maps. The likelihood 
function £(£) is concave, and in the present case it is 
defined on the convex set of CP maps. Its maximum is 
achieved by a single CP map if the data sample is suf- 
ficiently large, and the set of measurements is a quorum 
& 

The constraints to be imposed in the maximisation 
problem are the complete positivity and the trace- 
preserving property of the map £. A trace-preserving 
CP map is a linear map from operators in Hilbert space 
TL to operators in JC which can be written in the Kraus 
form fifll 



1 



(3) 
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Rjfy, F l (x l )=Aj(xi)A l {x l ) 



(14) 



where 



E44 = i«. 



(4) 



Let dim(7Y) = N and dim(/C) = M, and consider an 
orthonormal basis {U;} for the space of linear operators 
on Ti, namely 



and for any operator O 



N 2 



o = J2^o]Vi 



(5) 



(6) 



Upon defining the operator |ll| , [l2| 

N 2 



(7) 



i=l 



where * denotes complex conjugation, one can write for 
linearity 

£{p)^ r Tr H [{t K ®p T )S]. / (8) 
where T is the transposition. Notice that JljJ 

N 2 



X)Vi®V? = |*><*l , 



(9) 



where |\&) is given by the (unnormalized) maximally en- 
tangled state 



Hence, one has also 

£ = £®1(|¥)(*|). 



(10) 



(11) 



Eqs. (^) and (§1) establish an isomorphism between linear 
maps from 7i to K, and linear operators on the tensor- 
product space K £g> Ti. Complete positivity and trace- 
preserving property of £ imply 118] 



S>0 and Tric[S] = l n 



(12) 



For the construction of the likelihood function £(£) 
the condition S > is crucial JtJ. Actually, it allows to 
write 



S = C^C 



(13) 



where C is an upper triangular matrix, with positive di- 
agonal elements |l9| ]. Similarly, one has for the density 
matrices pj and the POVM's F t (xi) 



From Eqs. (g), (|13|) and (|14j), the likelihood functional 
in Eq. (0) rewrites 



A' 



C{£) ee C(C) =J2 lo S Tr i cic { R l R i 9A](xi)Ai(xi))] 
i=i 

K NM 2 

= E lo S X \{{n\C{R\ ® A\{ Xl ))\m))\ , (15) 

Z=l n,m— 1 

where denotes an orthonormal basis for Ti (£> /C. 

On one hand, the parameterisation in Eq. (]l^) implic- 
itly constrains the complete positivity of the map £ . On 
the other, the argument of the logarithm is explicitly pos- 
itive, thus assuring the stability of numerical methods to 
evaluate C{C). 

The trace-preserving condition is given in terms of the 
matrix S by Tr^S"] = 1^. However, the constraint 
Tt[S] = N which follows from Tr K [S] = ti H isolates 
a closed convex subset of the set of positive matrices. 
Hence, the maximum of the concave likelihood functional 
still remains unique under this looser constraint, and one 
can check a posteriori that the condition Tr^S] = 1-^ is 
fulfilled. Using the method of Lagrange multipliers, then 
one maximises the effective functional 



C{C) =C(C) -pTr[tfC] 



(16) 



where £(C) is given in Eq. (|T^) , and the value of the mul- 
tiplier p can be obtained as follows. Writing S in terms 
of its eigenvectors as S — J2i s i \ s i})(( s i\ > the maximum 
likelihood condition dC(C)/dsi = implies 

^■&[0>f®F,(*i))*il*i»«*i|] 



Tr[(pf ® F( Xl ))S] 



= /iTr[ Si | Si ))(( Si |] . (17) 
Multiplying by s.j and summing over i gives p = K/N . 



III. CHARACTERISING THE UNIVERSAL 
CLONING MACHINE 

We consider now the problem of estimating the CP 
map pertaining to the l-to-2 universal covariant cloning 
machine. The map is given by n 



£{P) = 3 s 2 {p<- 



1) S 2 , 



(18) 



where S2 is the projection operator on the symmetric sub- 
space, which is spanned by the set of vectors { | Si)(si \ ,i = 
0^2}, with | s ) = | 00), |s x ) = 1/V2(|01) + | 01)), and 
| s 2 ) = | 11), where {| 0), 1 1)} is a basis for each spin 1/2 
system. Using the lexicographic ordering for the basis of 
the tensor-product Hilbert space EG] one has 



2 



S2 



/ 1 \ 

1/2 1/2 

1/2 1/2 

V 1 / 



(19) 



In the following we label with A and B,C the Hilbert 
spaces supporting the input state and the two output 
copies, respectively. One can apply Eq. (^) to obtain the 
corresponding matrix S. Upon using the operator basis 
Vi = -^Vii with (To = 1, arid <7j (i = 1, 2, 3) denoting the 
customary Pauli matrices, one has 



5 = ^E[^) sc K B ®ac)( S2 ) sc ]^a; 



3 2 

= ~ [1a® (s 2 ) BC ] 

x [1a®( S2 ) SC ] 
2 

= - [1 A <g> (s 2 ) BC ] [\^)abab(^\ ® lo] 



L \i=0 



x [1a®(s 2 ) 



SCI 



(20) 



In the lexicographically ordered basis |2C 
He the matrix 5 writes 



of U A ®U B 



/4000022 
1 1 



S = 



(21) 



2 

1 1 2 
00000000 
00000000 
2 1 1 
2 1 1 
\0 2 2 4/ 



Now we want to apply the method presented in the previ- 
ous section to reconstruct the matrix S. We use random 
pure states at the input of the cloning machine 

- cos(0,/2)|O)a + e**» sin(^/2)|l) A . (22) 

In matrix notation pi = \ipi){ipi\ we write 



pi = -{t A + a A -ni) , 



(23) 



with ni = (sin 9i cos <pi , sin Qi sin </>; , cos 8i ) . On the two 
output clones we perform independent projective mea- 
surements along 
random directions rj = (sin a/ cos/3;, sin a; sin/3;, cos a;) 
and ti = (sin 7; cos 5;, sin 7; sin J;, cos 7;), then using the 
projector 

Fi{auh) = ~(1 B + a0 B ■ n) ® (l c + ha ■ U) , (24) 

where ai and bi denote the outcomes one has indeed ob- 
tained (the possible values are ±1). Defining 



at - 1 



a l = -+n- 4 



11 



at - 1 



7i = Y +7r " 4 



(25) 
(26) 



the Cholevsky decomposition for pf and i*/(<x/, 6;) writes 
as in Eqs. (Q), with 



/ cos(0,/2) e'>« sin(^/2) 
' 



(27) 



and 



A(at,bt) 



cos at e I/3 ' sin a; 











-id, 



cos 7j e ' sin 7/ 




(28) 



It follows that the matrix Qi(ai,bi) = Rj <X> Aj(ai,bt) 
that multiplies C in the likelihood function (|l^) is lower 
triangular with just the first column different from zero. 
One has explicitly 



[Qi(ai 
[Qi(ai 
[Qi(ai 
[Qi{ai 
[Qi(ai 
[Qi(ai 
[Qi{ai 
[Qi{ai 



M]n = cos(6»//2)cosa ; cos 7/ 
bi)]zi — e lSl cos(6*//2) cosq; sin 7; 
Mbi = e lA cos(6»//2)sina; cos 7; 
bt)] 41 = e l{l3l+Sl) cos(0i/2)sina; sin 7; 

i{0i/2) COS 5; COS 7; 



&Z)]51 — - ^^(^/WU^WJ ,( 

= e^'-^ sin(0i/2)cosaj sin 7/ 



6j)] 71 = e l(/3 '"*' ) sin(6' / /2)sin5 / cos7; 
6,)] 81 = e l(/3,+,5! -^ ) sin(6' / /2)sin5 ; sin 7, 



(29) 



We have now all the ingredients to construct the likeli- 
hood function in Eq. (|l5|), which will be a function of 
the 64 real parameters that specify the triangular matrix 
C. The problem of the maximisation of C(C) enters the 
realm of programming and numerical algebra optimisa- 
tion, where various techniques are known pl[ . 

In the following we show the results of a simulation 
obtained by applying the method of downhill simplex 
pl| , p2"t to find the maximum of the likelihood functional. 
This method is robust and efficient in case of a relatively 
small number of parameters. It has been reliably used 
in the reconstruction of the density matrix of radiation 
field and spin systems || , and in the characterisation of 
quantum communication channels for qubits H. 

The results are shown in Fig. [|. Pure states at the 
input of the cloning machine have been used, together 
with projective measurements over the two clones at the 
output. In both cases we adopted a uniform distribution 
on the Bloch sphere. The Monte Carlo method has been 
used to generate K — 10000 data, by using the theoreti- 
cal probability 



2 

p(ai,bt) = -Tr[s2{pi®l)s 2 Fi(ai,bi)} . 



(30) 
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6 S (i, j) 




FIG. 1. Maximum-likelihood reconstruction of the CP map 
of the l-to-2 universal covariant cloning. The picture repre- 
sents the values of the (real part of the) elements of the matrix 
S. Random pure states at the input of the cloning machine, 
and projective measurements along random directions on the 
two clones at the output have been used, with K = 10000 
couple of measurements. The statistical error in the recon- 
struction is of the order 10 -2 . The results compare very well 
with the theoretical values of Eq. 

A lengthy but straightforward calculation gives 

p(ai,bi) = - + 7(a/cosa ; + b t cos -fi ) cos 0j 
4 o 

+ \[ai sin ai cos(/3j - <j>i) + k sin 7; cos(<5; - <j>i)] sin 61 
6 

+— a; h [cosev/ cos 7; + sin a; sin 7; cos(5; — 0i)] . (31) 




the number K of simulated data. In accordance with 
the central limit theorem we find the asymptotic inverse- 
square-root dependence on K. 



IV. CONCLUSIONS 

In conclusion, we have applied a general method to 
reconstruct experimentally the completely positive map 
describing a physical device to the universal covariant 
cloning machine. The method, based on the maxi- 
mum likelihood principle, involves the maximisation of 
a functional, which depends on the results of quan- 
tum measurements performed on the clones at the out- 
put. The maximisation has to be made over all possible 
trace-preserving completely positive maps. A suitable 
parametrisation is allowed by the isomorphism between 
linear map from Hilbert spaces H. to K, and linear oper- 
ators in 7i (g) JC, along with the Cholesky decomposition 
of positive matrices. The numerical results we showed 
here has been obtained by applying the method of the 
downhill simplex to search the maximum of the likelihood 
functional. In our example, a good characterisation of 
the l-to-2 universal cloning machine has been achieved, 
with a number of simulated data as low as 10 4 . This 
is relevant, because some experiments are now feasible, 
but with low data rate or short stability time. In accor- 
dance with the central limit theorem, the statistical er- 
ror of the characterisation shows the inverse-square-root 
asymptotic dependence on the number of data. 

The method is very general, can be implemented im- 
mediately in the lab, and can be adopted in many fields as 
quantum optics, spins, optical lattices, atoms, ion trap, 
etc. 
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FIG. 2. Average statistical error in the characterisation of 
the universal covariant cloning machine versus number of data 
K. The error affects the value of the reconstructed elements 
of the matrix S that is univocally related to the CP map of the 
cloning machine. The dotted line represents the asymptotic 
dependence on the inverse square root of K, in accordance 
with the central limit theorem. 

In Fig. H we reported the value of the average statis- 
tical error that affects the matrix elements of S versus 



[1] W. K. Wootters and W. H. Zurek, Nature 299, 802 
(1982); H. P. Yuen, Phys. Lett. A 113, 405 (1986). 

[2] V. Buzek and M. Hillery, Phys. Rev. A 54, 1844 (1996); 
N. Gisin and S. Massar, Phys. Rev. Lett. 79, 2153 (1997). 

[3] D. Brufi, D. P. DiVincenzo, A. Ekert, C. A. Fuchs, 
C. Macchiavello and J. A. Smolin, Phys. Rev. A 57, 2368 



4 



(1998); D. Brufi, A. Ekert and C. Macchiavello, Phys. 

Rev. Lett. 81, 2598 (1998). 
[4] R. F. Werner, Phys. Rev. A 58, 1827 (1998). 
[5] N. Gisin and S. Massar, Phys. Rev. Lett. 79, 2153 (1997); 

N. Gisin and B. Huttner, Phys. Lett. A228, 13 (1997). 
[6] G. M. D'Ariano, C. Macchiavello, and M. F. Sacchi, 



quant-ph/0007062;; G. M. D'Arian o, F. De Martini, and 
M. F. Sacchi, |quant-ph/0012025| (to appear on Phys 



Rev. Lett.): G. M. D'Ariano and M. F. Sacchi, |quant- 



ph/000908q . 
[7] M. F. Sacchi, |quant-ph/0009104 . 

[8] S. L. Braunstein, A. S. Lane, and C. M. Caves, Phys. 

Rev. Lett. 69, 2153 (1992). 
[9] K. Banaszek, G. M. D'Ariano, M. G. A. Paris, and M. 

F. Sacchi, Phys. Rev. A 61, 10304(R) (2000). 
[10] G. M. D'Ariano, M. G. A. Paris, and M. F. Sacchi, Phys. 

Rev. A 62, 023815 (200 0); G. M. D'Arian o, M. G. A. 



Paris, and M. F. Sacchi, quant-ph/0009081 



[11] J. de Pillis, Pacific J. of Math. 23, 129 (1967). 

[12] A. Jamiolkowski, Rep. Math. Phys. 3, 275 (1972). 

[13] S. Yu, Phys. Rev. A 62, 024302 (2000). 

[14] J. I. Cirac, W . Diir, B. Kraus, and M. Lewenstein, |quant- 



ph/0007057. 



[15] For the concept of quoru m, see G. M. D'Aria no, L. Mac- 
cone, and M. G. A. Paris, |quant-ph/0006006] ; Phys. Lett. 
A 276, 25 (2000). 

[16] See, for example, M. A. Nielsen and C. M. Caves, Phys. 
Rev. A 55, 2547 (1997). 

[17] G. M. D'Ariano, P. Lo Presti, and M. F. Sacchi, Phys. 
Lett. A 272, 32 (2000). 

[18] For a positive — but not completely positive — map the 
condition of positivity of the matrix S is relaxed, by only 
requiring positivity for tensor product of vectors, namely 
(<f}\ ® (ip\S\<j>) <g> |V) > (see Ref. 

[19] Such decomposition is referred to as Cholesky decom- 
position, and is commonly used in linear programming 
(see, e.g., Ref. |2lJ). Moreover, if the matrix S is strictly 
positive the decomposition is unique. 

[20] We mean that \i) ® \j) precedes \k) <g> \l) if and only if 
either i < k or i — k and j < I. 

[21] See, for example, P. G. Ciarlet, Introduction to numerical 
linear algebra and optimisation, Cambridge Univ. Press, 
Cambridge, 1989. 

[22] W. H. Press, S. A. Teukolsky, W. T. Vetterling, B. P. 
Flannery, Numerical Recipes in Fortran: The Art of Sci- 
entific Computing, Cambridge Univ. Press, Cambridge, 
1992. 



5 



